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ABSTRACT 



Context. The cluster Abell 2163 is a merging system of several subclusters with complex dynamics. It presents exceptional X-rays 
properties (high temperature and luminosity), suggesting that it is a very massive cluster Recent 2D analysis of the gas distribution 
has revealed a complex and multiphase structure. 

Aims. This paper presents a wide-field weak tensing study of the dark matter distribution in the cluster in order to provide an alternative 
vision of the merging status of the cluster. The 2D mass distribution is built and compared to the galaxies and gas distributions. 
Methods. A Bayesian method, implemented in the Im2shape software, was used to fit the shape parameters of the faint background 
galaxies and to connect for PSF smearing. A careful color selection on the background galaxies was applied to retrieve the weak 
lensing signal. Shear signal was measured out to more than 2 Mpc (^ 12' from the center). The radial shear profile was fit with 
different parametric mass profiles. The 2D mass map is built from the shear distribution and used to identify the different mass 
components. 

Results. The 2D mass map agrees with the galaxy distribution, while the total mass inferred from weak lensing shows a strong 
discrepancy to the X-ray deduced mass. Regardless of the method used, the virial mass M200 falls in the range 8 to 14 10'"* /lygMo 
inside the virial radius (~ 2.0 h^g Mpc), a value that is two times less than the mass deduced from X-rays. The central mass clump 
appears bimodal in the dark matter distribution, with a mass ratio ~ 3:1 between the two components. The infalling clump A2163-B 
is detected in weak lensing as an independent entity. All these results are interpreted in the context of a multiple merger seen less than 
IGyr after the main crossover 

Key words. Gravitational lensing - Dark matter - Galaxies: clusters: general - Galaxies: clusters: individual (A2163) 



1. Introduction 

Clusters of galaxies are powerful probes for cosmology because 
they belong to the high-mass end of the halo mass function of 
collapsed structures dVoit 2005) . The mass distribution within 
these structures is representative of the nonlinear development 
of the accretion processes during the cos mological evolution 
JPress & Schechtei<1974HFrenket al.ll990l) . But it is still a chal- 
lenge to securely relate observational quantities like the X-ray 
luminosity and temperature, or the optical number counts and 
velocity dispersion, to constitutive prope rties such as the gravi- 
tational potential or the mass distribution (iReiprich & Bohringeij 
[2002; Vikhlinin et al. 200^, ICarlberg et all 119961) . Several in- 
ternal processes also play major roles throughout cosmologi- 
cal times. They can distort the simple scaling relations between 
physical quantities initially set by the hierarchical clustering pro- 
cess of structures formation. 

The rich cluster of galaxies Abell 2163 is an interesting clus- 
ter that displays several paradoxical properties. First, it is one of 
the richest Abell clusters, at a redshift z - 0.201, and among the 
most luminous o nes in X-ray s (Ly = 6.0 x 10"*^ h^^ erg s"' in the 
2-10 keV band. lElbaz et al.1 ( Il995h ). Extensive analysis of the 



global properties of this exceptional cluster have soon been con- 
ducted after the initial detection of a very high X-ray temperature 
(Tx - 14.6 ± 0.9 keV, measured with GINGA by Arnaud et aLJ 
(Il992h ). They have rapidly pointed evidence of a nonisothermal 
gas distribution with a s trong temperature drop in the outskirts 
dMarkevitch et al.lll994tlMarkevitchlll996h and inhomogeneities 
of the gas temperatu re in the ce nter ("Ma rkevitch & VikhlininI 
2001) . More recently, iFeretti et al. (2001 ) have detected a very 
extended and powerful radio halo, which they interpret as the 
tracer of a very hot nonthermal phase in the ICM. Moreover, ar- 
eas with a flatter radio spectral index and higher energetic parti- 
cles are found to be coincident with the most disturbed regions of 
the X-ray distribution (Feretti et al. 2004i) . Another component 
was identified as a secondary X-ray source in the map, spectro- 
scopically confirmed to be much cooler than the main cluster 
(A2163-B). The strong tempe rature variations in the center dis- 
played by Chandra imaging dGovoni et al.l 120041) are also rep- 
resentative of ongoing merging processes; however, the X-ray 
properties in the external regions of the cluster are more regular, 
both in luminosity and temperature, and fall within the general 
category of massive clusters. 
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Recently, iMaurogordato et al] (l2008h have presented a very 
detailed optical analysis of this cluster that sheds new light on 
its merging status. Using both spectroscopy of a large sample of 
galaxy members and wide field multicolor imaging, they stud- 
ied the galaxy distribution and confirm that the main cluster is 
dynamically separated from A2163-B, although they belong to 
the same complex. The galaxy density distribution of the main 
cluster A2163-A is strongly elongated in the EW direction and 
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shows a bimodal distribution that depends on the luminosity 
range of the galaxies tracing the distribution. For the faintest 
ones, the distribution is split into two components, Al and A2, 
for which the former is centered on the brightest cluster galaxy 
BCGl while the latter is slightly shifted co mpared to the galaxy 
BCG2 (following the naming defined by iMaurogordato et all 
(l2008h V In addition, a strong velocity gradient in the main clump 
is related to the spatial galaxy distribution and is elongated in 
the NE/SW direction. All these properties point towards a post- 
merging phase where the main component has just undergone 
a recent merger along the elongation direction, nearly copla- 
nar with the plane of the sky, while A2163-B is infalling onto 
A2163-A. Finally, a very detailed analysis of the complex X- 
ray emission, based on high-resolution imaging with Chandra, 
coupled with spectro-imaging with XMM -Newton, has brough t 
new light on the complex gas distribution (iBourdin et al.ll201lh . 
The authors clearly identify a cold gas clump in the direction of 
A2 but not coincident with it. This clump may be seen cross- 
ing the main cluster's hot-gas component, a scenario with fea- 
tures resembling those detected in the so-called Bullet cluster" 
IE 0657-56 (Markevitch et al. 2002). 

To explore the physical state of the cluster and to better un- 
derstand the relationship between the distributions of the galax- 
ies, the hot ICM gas, and the dark matter, we propose in this pa- 
per to use the weak gravitational lensing effect, which is directly 
related to the dark matter content and its distribution. A2163 
has already been observed in weak lensi ng by several groups 
but without strong and conclusive results. ISquires et al.l (Il997h 
used a large and shallow image of the cluster and found a mod- 
est value for the velocity dispersion of about 700 km s"' (ad- 
justed with a singular isothermal potential), although the authors 
claim that the uncertainties are so great close to the center that 
they could accept higher values up to 1000 km s ~'. With deeper 
data fr om the VLT but in a smaller field of view, ICvpriano et al.l 
(l2004l) measured o-y - 1020 ± 150 km s"'. Again this is a low 
value compared to what could be expected from such a hot clus- 
ter and t o what is measured dyn amically from the galaxies. More 
recentlv. lRadovich et al.l ( l2008l) have claimed that they solve this 
discrepancy with their new measure of the weak shear profile 
based on deep and wide field images of the cluster. However, 
their fit is not conclusive, so we decided to provide a new analy- 
sis of the weak lensing signal specifically dedicated to the spatial 
distribution of the dark matter, using the same CFHT/MegaCam 
data. We took special care in selecting the background galax- 
ies and remove the cluster contamination close to the center as 
much as possible. Moreover, our goal is to compare the 2D mass 
map with the X-ray map in order to characterize the link be- 
tween the dark matter and the baryonic mass better in this merg- 
ing cluster Similar approaches have led to the spectacu lar results 
in the "Bullet cluster" 1E0657-56 ( JClowe et al.ll2006l) in which 
the separation between the coUisional baryonic matter and the 
coUisionless dark matter is o bvious and allows a ccurate identi- 
fication o f the merging stage (iBradac et al.ll2006l) . In parallel to 
this work. lOkabe et al.l (1201 ll) used Subaru weak lensing data to 
provide a similar 2D analysis. Fortunately, their results are glob- 
ally similar to ours, although some discrepancies remain, and 
these are investigated in the paper 

The paper is organized as follow: in section 2 we describe 
the data and the weak lensing analysis, insisting on the process 
selecting the background galaxies. In section 3, an analysis of 
the central arclets is proposed and some strong lensing proper- 
ties are derived for the cluster center Section 4 presents the 2D 
weak lensing reconstruction and the spatial distribution of the 
dark matter in the field of the cluster, while in section 5 we dis- 
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Fig. 1. Galaxy number counts in the field of Abell 2163. The 
black dots correspond to r' magnitudes. The diamonds are the 
r' number counts in the CFHTLS deep field Dl. In both fields, 
only galaxies with CLASS_STAR < 0.8 are kept. The excess 
at bright magnitudes in A2163 represents the cluster galaxies. 
The limiting magnitude is consistent with what is expected for 
compact objects measured at the Scr level. 



cuss our mass estimates. Finally we present some conclusions in 
section 6. Throughout the paper we use Hq = 70kms ' Mpc"', 
Qm = 0.3, Qa = 0.7. At the redshift of the cluster (z = 0.2), 1" 
corresponds to 3.3 kpc and 1' to 200 kpc. Magnitudes are given 
in the AB system. 



2. Image data and weak lensing analysis 

2.1. Data reduction and photometry 

Imaging data were obtained at the Canada-France-Hawaii 
Telescope with MegaCam, during the run 05AC12 (PI: H. 
Hoekstra) and were retrieved from the CFHT archives at 
CADCQ. These data were alread y preprocessed by the CFHT 
Elixir pipeline (JMagnier & Cuilla ndre 2004) and correspond to 
Master detrend images. Nineteen images of 600 seconds each 
were found: 15 in r' and 4 in g'. The total integration time is 2.5 
hours in r' and 40 minutes in g'. 

Astrometric solutions for each CCD and then for each im- 
age were computed using SCAMP (Berlin 2006), a tool devel- 
oped at Terapijfl Absolute coordinates of the stars detected in 
the field were selected from the USNO-Bl catalog. Internal ac- 
curacy for the calibration of the whole set of images reached 
0.06", while the external accuracy was on the order of 0.27" 
or equivalently 1.5 pixels. All the images were combined and 
geometrically rescaled using SWARP, with a fixed pixel size of 
0.186". Final images in g' and r' cover one square degree. 

Photometric internal calibration is provided during the ob- 
serving run and photometric zero points are measured and in- 
cluded in each image header But some images were not taken 
in photometric conditions. They were post-calibrated with refer- 
ence to other avail able data, cross-corre l ating our catalogs with 
those provided by iMaurogordato et aP (l2008l) : these were ob- 
tained with the ESO Wide Field Imager on the 2.2m telescope 



' This research used the facilities of the Canadian Astronomy Data 
Centre operated by the National Research Council of Canada with the 
support of the Canadian Space Agency. 
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Fig. 2. Histogram of the average size ( yob) of the selected stars 
before (hatched histogram) and after (filled histogram) the de- 
convolution process implemented with Im2shape. For each ob- 
ject, the local PSF is determined locally by the average of the 5 
closest stars. 



and the Bessel R and V filters, under photometric conditions. In 
addition, they were corrected from the galaxy extinction, with an 
average value E{B - V) - 0.41. Therefore, the magnitude dif- 
ference between both sets of data, taking the AB correction of 
the MegaCam filters into account, defines the correction to add 
to the magnitude zero points: -0.08 mag in r' and -0.30 mag 
in g' . Both corrections include the galactic extinction correction 
and are estimated with an accuracy of 0.05 mag. 

Photomet ric catalogs were buil t in both filters, using 
SExTRACTOR dBertin & Amoutslll996h . The r' image, which is 
the deepest one, was the reference image for detection, and the 
g' - r' color index was computed inside the limiting isophote 
defined in the r' image for each object. Completeness magni- 
tudes were estimated from number counts to r'^.^^^ - 24.8 and 
s'tim ~ -^^-^ (Figure[T]i- They are approximately 0.7 brighter than 
the 5a- limiting magnitudes computed inside an aperture with a 
diameter 1 .45 times the full width at half maximum (FWHM) of 
the seeing disk, and defined in the MegaCam exposure time cal- 
culator (DIETIfl). This last magnitude is best suited to point-like 
objects, and that the 0.7 mag difference reflects that most of the 
faint objects are compact but extended ones. In the final images 
used in this paper, the measured average seeing is 0.75" in r' and 
0.67" in g'. 



2.2. PSF correction and galaxy shape measurements 

Galaxies are first separated from stars with the help of the 
magnitude-peak surface-b rightness diagram, wh ich clearly iden- 
tifies the stellar objects dBardeau et al.l l2005h . Then for each 
galaxy, the local PSF is measured by averaging the shape param- 
eters of the five closest stars. Intrinsic galaxy shape parameters 
are finally recovered using the Im2shape software developed by 
[Bridle et al. (2001). This software performs an analytical decon- 
volution of the galaxy images, approximated by a single ellip- 
tical Gaussian, and returns the intrinsic shape parameters with 
their uncertainties, estimated by a Bayesian analysis of the im- 
age residuals. It has been successfully applied on a sample of 
clusters of galaxies by .Bardeau et al. (2007i) . who developed a 
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Fig. 3. Color-magnitude plot of the galaxies in the field of Abell 
2163, located within a 10' radius from the central galaxy BCGl. 
The straight line corresponds to the cluster red sequence (RS), 
and the dashed lines define the window for selecting "cluster 
members". Galaxies located less than 200" from the cluster cen- 
ter and within the selection window are plotted as large dots. 



methodology for the shear analysis that we follow fairly closely. 
To test the validity of the analysis chain, we checked that the re- 
sult of the deconvolution process applied to the stellar field itself 
shrinks the size of the deconvolved stars to less than 0.1", with a 
flat distribution of their orientation (Fig.|2}. For the galaxies, the 
histogram extends to a larger size, up to several arc seconds. The 
whole pipeline was also applied on th e simulated sheared im ages 
developed by the STEP consort ium (iHevmans et al.ll2006h . The 
results are presented in detail in iFoex et al.l (l2012l) . In short, us- 
ing simulated sheared images of faint galaxies, the average mea- 
sured shear is compared to the true shear input in the simulations. 
This is done on a sequence of simulations with increasing shear 
value and several sets of realistic PSFs. The average difference 
bewteen the shears is fit with a linear function of the true shear. 
With our pipeline we systematically underestimate the true shear 
by 10%, but with no systematic offset in the measured shear. 
Consequently, we corrected the measures of the shear by apply- 
ing a factor 1 . 1 to the value of average tangential ellipticity of 
the galaxies. This correction affects the total mass determination 
by ~ 15% upward. 

2.3. Selection of background galaxies 

An important step in building a photometric catalog of objects 
for weak lensing is to clean it from cluster and foreground con- 
tamination as much as possible. The cluster contamination is the 
major source of errors in the mass reconstruction for low red- 
shift clusters: it is most prominent in the center so it distorts the 
shear profile and attenuates its central value, leading to an un- 
derestimation of the central mass density. In practice, we first 
selected the galaxies within a magnitude range large enough to 
have a high number of objects but limited in the bright end to re- 
move a large fraction of cluster members. The lower limit of the 
magnitude cut was chosen as r' - 22, a magnitude fainter than 
that of a 0.1 L* galaxy at the cluster redshift. The upper limit was 
fixed to r' - 25.5, i.e. 0.7 magnitude fainter than the complete- 
ness limit. This limit corresponds to the 5cr limiting magnitude 
of point sources estimated with DIET. We added a selection cri- 
rleBjgJdibaaediaBaMthhtiiiai or index. This kind of selection by the 



G. Soucail: Dark matter distribution in thie merging cluster Abell 2163 



uuuu 




' 


' ' 






z (cluster) 


- 


8000 


- 




p 




- 


6000 


- 


^ 








- 


4000 


- 






-1 


- 


2000 


- 




1 : 









^^^^_^^^^ - 



12 3 4 

photometric redshift 

Fig. 4. Photometric redshift distribution of the "reference cat- 
alog" built from the CFHTLS-Dl data. The photometric se- 
lection applied to the weak lensing catalog has been included 
(see text for details) as well as a photometric redshift cut with 
Zzphot > 0.35. 



that have been publicly available after the T0004 release. They 
were carefully calibrated and validated with spectroscopic sam- 
ples (Coupon et al. 2009). We applied to this catalog the same 
selection criteria as we applied to the present catalog (magnitude 
cut 22 < r' < 25.5 and color selection with g' -r' outside the red 
sequence window), and we only selected those objects having a 
photometric redshift Icr error smaller than 0.15. We assumed 
that this sample is representative of the galaxy population on the 
line of sight of A2163. Looking at its redshift distribution, we 
measured a foreground contamination of 4%. Moreover, for all 
the galaxies with a photometric redshift Zphoi > Zduster + 0. 15, we 
computed the average geometric factor for weak lensing (Fig|4]i: 



Dr 



Df 



>= 0.72 + 0.10 



with an average redshift for the sources < Zphot >= l-l- Our 
weak lensing catalog includes objects fainter than the complete- 
ness limit (r' ~ 24.8) so a significant fraction of the galaxies is 
probably missing in the range [25 - 25.5]. But we checked that 
it does not significantly affect the redshift distribution, and the 
average factor < Dis IDqs > is not changed by more than a few 
%. 



galaxies color propertie s has already been tested in many stud- 
ies (Kneibetal. 2003c iBroadhurst et al.l 120051: [Limousin etall 
l2007h . It has proved to be quite powerful to remove galaxies 
with photometric properties compatible with cluster members. 
We built the color-magnitude diagram from the photometric cat- 
alog and identified the cluster red sequence (Fig. 13). The average 
color index of the bright ellipticals is g' - r' - 1.15 + 0.05 mag, 
in goo d agreement with expecte d values from galaxy evolution 
codes dBruzual & Charlotll2003l) . All galaxies located +0.15 and 
-0.20 from the red sequence line were stored in a catalog that 
we identify as the "cluster galaxy" catalog. Although it is not 
fully representative of the color diversity of cluster members, it 
includes early-type galaxies that dominate the density distribu- 
tion of the galaxies. Moreover, a fraction of bluer ones is also 
included because the lower color limit is slightly widened. The 
galaxies outside this window and within the magnitude range 
[22.0-25.5] are stored for the weak lensing analysis. The result- 
ing catalog contains 14550 galaxies within a square area of side 
37' centered on the galaxy BCGl. This sample is expected to 
be substantially free of cluster contamination and is dominated 
by background sources. With this selection, the mean density of 
sources reduces from 18 to 12 galaxies/arcmin^ in the central 
area of the field, with a residual slight excess compared to the 
most external regions of the field (no more than 10%). Thus this 
corresponds to an efficient cleaning of the cluster contamination. 
The quantitative mass reconstruction based on the measure 
of the shear signal also requires an absolute scaling by the lens- 
ing factor < Dis IDqs >, which depends on the redshift distribu- 
tion of lensed galaxies. In this ratio, D^s is the angular diameter 
distance between the lens an d the source while Dp s is between 
the source and the observer (ISchneider et al.l[l992l) . We do not 
have direct information on this redshift distribution in the present 
case, so we attempted to estimate it carefully. To do so, we built a 
photometric catalog considered as a reference catalog, using the 
T0004 release of the CFHTLS-Deep survey. Its main advantage 
is that the data have been collected with the same instrument and 
the same filter set. But they are much deeper than the present data 
(at least one magnitude deeper in g' and r') and multicolor ob- 
servations in five filters allow the determination of accurate pho- 
tometric redshifts. In practice we used the photometric red shifts 



3. A two-dimensional mass map 

To map the 2D mass distribution, we use the LensEnt2 method, 
an entr opy-regularized maxim um-likelihood technique devel- 
oped bv [Marshall et al.l (|2002|) . The mass map is built from the 
weak lensing catalog, using a smoothing scale characterized by 
the so-called intrinsic correlation function (ICF) of the model. 
A Gaussian ICF of 160" (~ 520 h^^ kpc at the cluster redshift) is 
a good compromise between details and smoothness of the mass 
map, given the number density of sources in the background cat- 
alog. The final mass map is scaled in units of mass surface den- 
sity {MqPc'^) and can be used directly to estimate the masses of 
the different mass components. The code also generates an error 
map that allows the signal-to-noise ratio map to be built and the 
significance of the detected peaks to be quantified. 

In the field of Abell 2163, the main peak is detected at more 
than 8cr and is centered close to the central galaxy BCGl. No 
strong lensing features or multiple images are identified around 
this galaxy, precluding any strong lensing model to explore the 
central mass distribution. Howe ver, a few single imag ed arclets 
were reported around the BCGl (iFort & Mellierll994l) . confirm- 
ing that a peak of dark matter coincides with the galaxy (Fig.|6]l. 

The mass distribution is also significantly elongated towards 
the second brightest cluster member BCG2 (Fig. |5) with a sec- 
ond mass peak detected at ^ 6o". A detailed comparison between 
this mass map and the galaxy density obtai ned with the galaxies 
selected within the "red sequence" (Fig. 6 in lMaurogordato et al.l 
(2008)) shows a striking concordance between the two distribu- 
tions. Even the much less significant peaks in the mass distribu- 
tion are coincident with other subclusters, namely the northern 
subcluster A2163-B, and other peripheral clumps C and D. The 
correlation between the dark matter distribution and the cluster 
galaxies density distribution clearly demonstrates that they both 
trace a noncollisional component in merg ing clusters. These 
maps are very differ ent from the X-ray map (iGovoni et al.ll2004t 
iBourdin et al.ll201 lb . which traces the gas distribution with the 
following characteristics: a single peaked distribution, with an 
E-W elongation and a peak center significa ntly off-center com- 
pared to BCGl . iMaurogordato et al] (l2008b propose in their pa- 
per a scenario in which we observe A2163 in a post-merger 
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Fig. 5. MegaCamr' image of the cluster Abell 2163, overlaid with the mass map reconstruction. The lowest contour corresponds to 
the 2cr level and the following ones are scaled with cr. The clump A2163-B is detected at the 3.5cr level. 



phase, less than 1 Gyr after the main collision. At this stage 
the gas morphology is expected to still be centrally concen- 
trated, and elongated along the merger axis with a possible dis- 
placement of the gas relative to the main dark matter clump to- 
wards the secondary clump. This scenario remains fully compat- 
ible with the present observation of the dark matter distribution 
traced by weak lensing and is verified by numerical simulations 
of cluste r mergers including bo th dark matter and intracluster gas 
physics dRoettigeret al.|[l997h. Moreover, it has b een recently 
confirmed by the analysis of iBourdin et al.l (1201 ll) . who show 
evidence of a cool gas "bullet" close to the cluster center that is 
clearly separated from its galaxy component and identified here 
as the second mass peak of the cluster. 

As already mentioned, the northern secondary peak is also 
detected well in the mass map, at the 4cr level. It corresponds 
closely to the clump A2163-B identified by Maurogordato et al. 
(l2008l) as an inf ailing group with a separated velocity struc- 
ture. However, there is some uncertainty in this mass distribu- 
tion, which seems strongly elongated in this area. It is difficult 
to know whether this is related to an artifact in the mass re- 
construction or to a realistic mass extension. The proximity of 
a bright star close to the mass clump (3' north of the clump, see 
Fig. Is) reduces the available area where background galaxies 



can be used for shear measurements. Moreover, iDietrich et al.l 
(1201 ll) show that peak offsets in weak lensing maps can occur up 
to 1' when the reconstruction is dominated by shape noise, i.e. 
the statistical noise due to the intrinsic ellipticity of the sources. 
This is indeed the case in this low S/N area of the mass map. 
The comparison w ith the galaxy density distribution (Fig. 5 in 
iMaurogordato et al. (2008) ) also shows similarities in this elon- 
gated distribution, so regardless of the exact centering of the 
weak lensing mass clump, it is closely associated with the X- 
ray peak and the group A2163-B infalling on the main cluster 
We therefore confirm that a significant amount of mass is associ- 
ated with A2163-B. We propose to quantify better it in the next 
section. 



4. Masses in Abell 2163 

4. 1 . Weak lensing masses 

There are several ways to estimate the mass from the shear mea- 
surement of the background galaxies. The simplest one is to fit 
the shear profile directly by some analytical profiles, like the one 
expected from a singular isothe rmal sphere (SIS) or an NFW dis- 
tribution (iNavarro et al.lll997l) . These mass profiles correspond 
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Fig. 6. Combined 3 -color image of the center of the cluster Abell 
2163 obtained with a combination of V, R and I images taken 
with FORSl at the VLT in exceptional seeing conditions of 
0.50". In this image, the identification of the most significant 
highly distorted features is reported. The central galaxy BCGl 
is located at: aoooo = 16''15"'49.0'; ^2000 = -06°08'42.8". 
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Fig. 7. Measured shear profile, averaged in sliding windows (cir- 
cular annularrings of 1 10" each). Data with 100" < r < 250" or 
375" < r < 450" have been removed from the fit because radial 
signal may be perturbed by the additional mass clumps A2 or B. 
The best fit with an SIS mass profile and an NFW one are plot- 
ted. The fit is limited to 700" around the cluster center because 
the signal is not significant enough at a greater distance. 



respectively to the density profiles 
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where pc is the critical density at the cluster redshift. 6c is 
related to the concentration parameter c,.,, by 
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Fig. 8. Distribution of the fitted parameters of the NFW profile in 
the M200 - c,.,> plane, scaled by their likelihood. The coiTelation 
between both parameters found from numerical simulations of 
dark matter halos is plotted as a dashed line. 



and the virial radius r2oo is defined by r2oo = chirrs- The virial 
mass M200 is the mass enclosed inside the virial radius: 
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For all the fits, we assume spherical symmetry. However, we 
have shown in the previous section that subclustering is quite 
significant in the mass distribution, so we decided to remove 
from the fit the data points located inside two concentric cir- 
cular annuli corresponding to A2163-A2 (100" < r < 250") 
and A2163-B (375" < r < 450"). We also shifted the center 
of the lens 1.3' west from BCGl and fixed it between the two 
main mass peaks, to enclose the whole mass distribution better. 
In practice, this shift does not have a strong impact on the shear 
profile at large radius. The shear profile is shown in Fig. with 
the result of the best fits by the two mass profiles. In both cases, 
we limited the extent of the fit to 700" because we consider that 
the signal is not reliable at a larger distance. 

Instead of binning and averaging the shapes of the distorted 
galaxies in circular annuli, we also implem ented a global ap- 
proach. We us ed the McAdam software fiMarshaU et al.l I2OO2I 
lMarshallll2006h . which is based on a Bayesian analysis of the 
shear distribution, approximated by a given parametric solution. 
The output of McAdam is a probability distribution of the fitted 
parameters, obtained using a maximum-likelihood estimator and 
an MCMC iterative minimization. McAdam works directly on the 
PSF-corrected faint galaxies catalog. It takes all the galaxies into 
account at their position with their shape parameters, without 
any ring averaging. Again the objects close to A2163-B are re- 
moved from the working sample, and the mass center is shifted 
1.3' west of BCGl, with the goal of estimating the global com- 
ponent of A2163-A. 

Thanks to this global fit procedure, the shear field was fit 
with the SIS and NFW distributions, assuming spherical sym- 
metry in both cases. A prior on the concentration parameter c 
is included (1.5 < c < 10), following the results of N-body 
simulations o f cosmological structure formation at the galaxy 
cluster scale (iBuUocket al.ll200U iHennawi et"aDl2007h . Fitting 
the shear profile with the SIS mass profile yields a velocity dis- 
persion cTsiieai or, equivalently, a value for the Einstein radius 6^ 
scaled by the ratio Dls/^^os averaged over all the sources. For 
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Fig. 9. Integrated mass profiles obtained with the best-fit param- 
eters of the SIS and NFW density profiles. The values of the 
projected mass directly measured on the 2D mass map and in- 
tegrated up to the radii 1 Mpc and 2 Mpc (^ r2oo) are plot- 
ted for comparison. They correspond to projected masses inside 
cylinders so they overestimate the mass by ~ 25%. The total 
mass also includes the contribution of A2163-B with a mass of 
~ 2.7 10'^ /z^Q Mq. The vertical dashed line corresponds to the 
position of the virial radius raoo- 



the NFW mass profile, there is a well known degeneracy be- 
tween M200, the mass enclosed within the virial radius r2oo, and 
the concentration parameter c,,„ (Fig.|8]l. But the average value of 
Cvir resulting from our fit (c,,,> = 4.6 + 1 .6) is consistent with the 
values pre dicted by numerical simulations of dark matter halo 
properties (iBuUock et al.ll200l1) . The results of the different fits 
are presented in TablelT] and the integrated mass profiles are dis- 
played in Fig|9] 

The masses of the different mass clumps identified in the 2D 
mass map (Fig.jSj are also measured through direct integration of 
the mass distribution. No special shape is assumed and the mass 
is integrated in hand-defined contours or within given apertures. 
In a first attempt, we integrated the mass of the main clump in 
concentric annuli up to 300"(- 1 Mpc) and centered between 
the two mass peaks Al and A2. This radius corresponds to the 
limit where A2163-B starts to be detected in the map, and it in- 
cludes both mass peaks (Al and A2). We find a projected mass 
of 7.3 10'"^ /i^Q Mq, which is the mass integrated within a cylin- 
der of radius 1 Mpc. This projected mass exceeds the mass inte- 
grated inside a sphere of the same radius by approximately 25% 
for a standard NFW profile. In comparison, the projected mass 
deduced from the previously fitted values of the NFW profile 
and enclosed in a 1 Mpc radius is 8.0 + 2.3 10'"* /i^q Mq, con- 
sistent with the direct value. We also measured the total mass 
inside a cylinder defined by the "virial radius" r2oo taken as 2.0 
Mpc. This mass ~ 1.7 10'^ /i^,} Mq includes the contributions 
of the two mass clumps, A2163-A (Al and A2) and A2I63-B. 
It reduces to M,o, ~ 1.4 10''' /i^,] Mq for the 3D mass inside 
a sphere of radius ^200 centered between Al and A2. This last 
value is considered as the total mass of the cluster within the 
virial radius (Table[T). 

Going further in the mass measures from the 2D mass map, 
we also attempted to estimate separately the masses of each of 
the two mass subclumps that we have identified in A2163-A (Al 
and A2). These are approximate measures because they are spa- 
tially limited and do not properly include the large-scale exten- 



sion of the cluster The measured masses are 7.1 10'"'/i.jQ M0 and 
2.5 lO'"* /i^Q Mq respectively, and their sum represents a major 
fraction of the total mass. Interestingly, we are able to measure 
a mass ratio ~ 3:1 between the 2 clumps, although this result 
must be taken with cautio n. This is not too dif ferent from the 
mass ratio 4: 1 estimated bv lBourdin et al.l (1201 ll) from interpre- 
tating of the merging scenario. Our ratio 3:1 also seems closer to 
the apparent ratio in the galaxy distribution and could be more 
representative of the total mass repartition. 

For the clump A2163-B, we integrated the mass within a 
hand-designed contour close to a circular annulus of 500 kpc, 
in order to avoid the contamination by the main cluster This 
gives a mass of ~ 2.7 lO'"* /j^q Mq, so sUghtly higher than the 
mass estimated bv lBourdin et al.l(l20I ih . However this is a rather 
uncertain measure, since the weak lensing signal and its distri- 
bution are not fully reliable. 

All these results are summarized in Table [T] It is interesting 
to note that the sum of the direct mass of the individual clumps 
(Al + A2 + B) represents ~ 75% of the total mass. The remain- 
ing mass represents the large-scale extension of the distribution 
and possibly the additional clumps contribution. Moreover, al- 
though there are large uncertainties in the direct approach, the 
results are consistent with those derived from the fit with usual 
mass distributions. 



4.2. Comparison with other lensing mass measurements 

It has long been pointed that, regardless of the method used to 
determine the mass of A2163, the lensing estimates are much 
lower t han the values expec ted for this very hot cluster For ex- 
ample, ISquires et aP (Il997h fitted their data by an SIS profile 
with a velocity dispersion ctsis - 740 km/s, and lCvpriano et all 
(2004) found cr = 1020 ± 150 km/s. Both values are compatible 
with the present measure, even if the y were obtained with fe wer 
galaxies: 261 bac kground galax i es for l Cvpriano et al1(l2004l) and 
736 galaxies for ISquires et al.l (1 19971). We are more s keptical 
of the fitted values proposed by iRadovich et al.l (l2008h derived 
from the same initial data as in the present work. With our se- 
lection process and the great care we took for the PSF correc- 
tion, we ended with a higher galaxy density (12 instead of 8 
galaxies/arcmin^), providing the largest sample of background 
sources used in the weak lensing mass determination in A2163 
14550 galaxies within 0.38 squ are degree). Moreover , looking 
into details in the shear profile in lRadovich et al.l(l2008l) . we note 
that their fit starts with data points located at a much larger ra- 
dial distance than ours. From our study, removing data points 
below r < 200" in the fit of the SIS profile inc reases the value 
of cr t o 1 100 km/s, close to the value found by IRadovich et al.l 
(I2OO8I) . The quality of the fit is very poor so we consider that 
the large-scale mass distribution in A2163 is too complex and 
shallow to be fitte d by a single SIS profile. More recently, 
lOkabe et al.l (l20Ilh have also performed a weak lensing analy- 
sis using Subaru imaging data. They confirm the bimodal struc- 
ture of the mass in the central parts of the cluster, although the 
mass ratio between Al and A2 is much higher (8:1 to 10:1), 
and the detection of the clump A2 is not as significant as in our 
own map. The 2D mass extractions from weak lensing recon- 
structions are still Hmited in their quantitative use and would 
strongly benefit from the h igher source density of HST imaging 
("Becker & Kravtsov|| 20I ih . The total mass estimates obtained 
by Ok abe et al. ( 201l1) are also higher than the present ones, at 
least by 50% when they use a multicomponent analysis and by 
a factor 2 for a single mass profile reconstruction. But in the last 
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Table 1. Best-fit values obtained from the fit of the SIS and NFW profiles. 



NFW 

^■200 (/i7o Mpc ) M200 (10'" /zto' Mo) rjoo (/i7o Mpc ) M500 (10'" fco' A^o) 



Radial bins 
Global fit 



3.36: 
4.6: 



0.09 
1.6 



1.97: 

1.79 i 



:0.2 

0.23 



10.7: 

8.0 d 



:3.1 

3.0 



1.27: 
1.18: 



0.12 
0.13 



7.1: 

5.7: 



2.1 
2.1 



eE (") 



SIS 

o-ios (km/s) 



Mm (r < rjoo) 



Radial bins 
Global fit 



13.7: 
13.0: 



3.6 
2.6 



810±110 
790 ± 80 



600 ± 220 
570 ±180 



2D mass map 



M^D (r < 1 Mpc) ^ 7.3 10" h-^^ M^ 

M2° (A2163 - Al) ^ 7.1 10'" /i,^' M^ 

Nf^ (A2163 - A2) ^ 2.5 10'" h-^, M^ 

M-° (A2163 - B) ^ 2.7 10'" h-^ Mq 

M^^ (r < r2oo) ~ 17.0 lO'" hj^ Mq <:==> M, 



"70 



14.0 10'" h-^ Mo 



The first line corresponds to the ID fit of the shear profile while the second line gives the result of the global fit with McAdaiti. The estimates based 
on the strong lensed features are indicated as a complement. The last lines correspond to the projected masses measured directly on the 2D mass 
map. They correspond to masses measured inside a cylinder of radius 1 Mpc and 2.0 Mpc (= r2oo), respectively, as well as the projected mass of 

the different mass components A2163-A1, A2163-A2 and A2163-B. 



case, their fit show high mass associated with a low value of the 
concentration parameter (c ~ 2.8). This combination of fitted 
parameters is probably representative of the mass-concentration 
degeneracy shown in Fig |8] 

In conclusion, there are still some discrepancies between 
different approaches used to infer quantitatively the total mass 
of this complex cluster, but the more recent results become 
more consistent with each other, provided similar criteria are 
taken into account to define a "total mass". Absolute values for 
the mass still suffer from uncertainties in the large-scale struc- 
tures' contributions and are also highly dependent on the intrin- 
sic noise owing to th e limited number of background sources: 
iHoekstra eTail (1201 Ih show that it can reach a 20 % to 25 % 
level for ground-based observations, whatever the methodology 
established to produce the results. 



4.3. Galaxies, gas and dark matter distributions: 
consequences on tlie merger scenario 

We focus in this section on a comparison between the masses 
derived from weak lensing, which trace the dark matter in the 
cluster, the masses derived from X-rays, and the spatial dis- 
tribution of the galaxies. The X-ray temperature is very high 
{Tx ~ 12 to 14 keV) regardless of the instrum ents used to 
measure it (lEIbaz et al.lll995l: iGovoni et aT1l2004 . The temper- 
ature distributi on shows strong varia t ions in the central part 
of the cluster dMarkevitch et al.|[T996t iMarkevitch & Vikhliiiinl 
l200lh and a significa nt nonthermal component at high energy 
(iRephaeli et al.ll2006h . Extrapolating the mass under these con- 
ditions is a risky exercise and must be taken with caution be- 
cause the cluster is obvi ously not in t h ermal equilibrium. To 
overcome this difliculty, iBourdin et al.l (1201 1|) hav e computed 
the Yx parameter defined by Kravts ov et al.l (|2006l) . Yx is con- 
sidered as one of the best proxies for the total cluster mass. 
They used different calibrations of the scali ng relation Yx - M500 
(lArnaudetaI.ll2010t IVikhUnin erai]|2009t) to give a total mass 
estimate that does not rely explicitly on the hydrostatic equilib- 
rium hypothesis. Regardless of the calibration used, they find a 
value Msoo - 1.8 to 2.0 10'^ Mq, i.e.. This converts to a mass 
M200 within the virial radius 50% higher for a typical NFW mass 



profile, so this mass is two to three times higher than the lens- 
ing value. It is not clear why there is still such a significant dis- 
crepancy between the weak lensing mass estimates and the ones 
derived from X-ray proxies. 

The spatial distributions of the three main components 
of the cluster are d isplayed in Fig [TO] the intrac luster gas 
(iBourdin et al.ll201 ll) and the dark matter (this paper) distribu- 
tions are overlaid on the galaxy cluster image. As already ob- 
served in several other merging clusters, these distributions show 
evidence that the baryonic mass (the gas) does not follow the to- 
tal mass (dominated by the dark matter). On the contrary, the 
galaxies distribution is spatially more coherent with the dark 
matter distribution A bi-dimen sional Kolmogorov-Smimov test 
dFasano & Franceschinil[l987l) could be applied to get a more 
quantitative answer but his needs further adjustments in the 
galaxies selection, which are beyond the scope of this paper. 
[Maurogordato et al] (l2008l) show morphological differences in 
the galaxy density distribution for different luminosity bins that 
require deeper insights to reconcile with the dark matter distribu- 
tion. Nevertheless, in contrast to the intracluster gas, both galax- 
ies and dark matter represent noncollisional mass components 
so we confirm that the total mass of the cluster is dominated by 
a collisionless form of dark matter This is at least the fourth ex- 
ample of a cluster merger with such proper ties, after the initial 
work on the "Bullet cluster" IE 0657-56 (iBradac et al.ll2006h 
followed by similar trends observed in Abell 520 ('Ma hdavi et al.l 
I2OO7I) and MACS J0025.4-1222 (Bradac et al. 2008) . 

Abell 2163 is the cluster with the strongest discrepancy be- 
tween X-rays and lensing masses. It is hard to reconcile both 
values, even when knowing the different limitations discussed 
above, and we still do not have a clear vision of the physical 
origin of the discrepancy. Because of the merging process, the 
gas distribution could be far from isothermal on large scales, 
and this might induce a strong uncerta inty, which depen ds on 
the exact phase of the cluster merger (iPoole et al.l l2007h . We 
suspect that all absolute mass determinations with X-ray prox- 
ies in a s trong merger like A2163 must be taken with caution 
dOkabe et al. 2010). and the gas mass must be revisited in the 
framework of a colliding cluster We also expect better lensing 
mass estimates using deep HST data. However, apart from a dif- 
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Fig. 10. Multi-color image of the center of the cluster Abell 2163 with overlays of the dark matter distribution (green contours from 
the lensing map) and of the gas distribution (red contours for the Chandra X-rays image, courtesy H. Bourdin). The shift between 
the two mass clumps is 3', while the distance between the X-ray maximum and the eastern clump A2163-A1 is 1.2'. The total size 
is lO'x 10'. 



ference in the absolute mass scaling, the relative distributions of 
the three components are consistent with each other in a general 
scenario of a cluster major merger. 

5. Conclusion 

We have presented in this paper a detailed dark matter map of the 
cluster Abell 2163, using the weak lensing effect on background 
sources. We provide global masses, as well as the bi-dimensional 
mass distribution and mass estimates for each individual sub- 
structure identified in the map: 

- Thanks to a detailed mass reconstruction, we were able to 
separate the central mass clump into two components with 
a distribution in agreement with the clus ter galaxies den- 
sity distribution jMaurogordato et al1l2008b . The mass of the 
main clump is about 1.0 lO'^ /z^g Mq and is split in two 
clumps with a mass ration ~ 3:1. It is consistent with pre- 
vious weak lensing estimates, but does not follow the ex- 
pected correlations with the X-ray properties of the gas. Such 



a mass discrepancy is a strong argument in favor of the merg- 
ing scenario in the cluster. The second mass clump A2163- 
A2 is shifted comp ared to the cold gas core detected by 
iBourdin et al.l (1201 ll) . and this shift can be due to the trail- 
ing of the gas compared to the noncollisional dark matter 
during the crossing of the main clump. 
We detect at more than 4cr the dark matter associated 
with the clump 2163-B, identified spectroscopically by 
iMaurogordato et aP (l2008l) . This clump has a mass value 
(^ 2.7 lO'"* /z^g Mq) compatible with a small cluster. This 
is also coherent wit h the X-ray temperatu re on the order of 
4 keV identified by Bourdin et al. ('20111). As already sug- 
gested by Maurogordato et al. (2008), this clump has not yet 
undergone its main interaction with the central cluster and is 
still infalling towards the center. 

Concerning the total mass of the whole complex, our mass 
estimates do not exceed 1.5 10'^ /i^g Mq. This value is well 
below the expected mass given by the coiTelations with the 
X-ray properties of the cluster. We confirm that the scaling 
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laws that use X-ray proxies to infer the mass of clusters must 
be taken with caution for structures that are far from the hy- 
drostatic equilibrium. 

The results of our study are a new step forward in the pro- 
cess of building a coherent scenario for the multiple mergers 
in the cluster Abell 2163. This thorough analysis of the weak 
lensing mass reconstruction was possible thanks to the combi- 
nation of both the spatial extent and the depth of the data. It 
confirms the bimodal structure of the matter distribution in the 
cluster center and the physical separation between the baryonic 
mass and the noncollisional mass components. Our results also 
show that there exists many similarities between the gas and dark 
matter distributions in Abell 2163 and in the "Bullet cluster" 
1E0657-56. Both cases display a gas core spatially separated 
from its dark matter component, even if the mass ratio between 
the clumps, the 3D configuration of the merger and the age of 
the merger differ. But they offer different merging configura- 
tions that can feed the comparisons with numerical simulations 
dPoole et al.ll2007h . Moreover, Abell 2163 presents an additional 
complexity with the identification of several other mass clumps 
involved in the merger but at an earlier stage or in the premerging 
phase. In conclusion, Abell 2163 has long been considered as an 
exceptional cluster because of its high gas temperature. It is now 
more interesting to explore the richness and the complexity of 
the physical processes occurring within it. 
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